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O ABSTRACT 
O 



O 



We present Ha rotation curves of two low luminosity field galaxies with 

-0.8 , j it* _ 1/1 <7+0. 3 



r-band absolute magnitudes M r = -13.9Iq:5 and M r = ~^.7Z 02 (for h 



// u 100 km s 1 Mpc = 0.7; the large error bars reflect distance uncertainties). 
Most previously studied galaxies in this luminosity range are members of groups 
defined by brighter galaxies, but these two systems, selected from Blanton et al.'s 
(2004) sample of low luminosity galaxies in the Sloan Digital Sky Survey (SDSS), 
appear to have no bright companions. The measured rotation speeds at the outer 
extent of the Ha rotation curves, 34.8±3.8 kms -1 and 30.9±7.2 kms -1 , are larger 
than the values of I6.OI5 4 kms -1 and 20.9^52 kms -1 predicted by extrapolating 
the inverse Tully-Fisher relation of luminous SDSS galaxies to these faint lumi- 
nosities. However, a previous HI measurement of the first galaxy shows that it 
has a gas mass similar to its stellar mass, and the total baryonic mass is consistent 
with that predicted by McGaugh et al.'s (2000) "baryonic Tully-Fisher relation." 
We find r-band dynamical mass-to-light ratios within the radii of the last Ha 
data points (about 1.8 disk scale lengths in each case) of 12.6^1 M & /L Q and 
4.8l|i Mq/Lq, much higher than the values ~ 1 M /L expected for the stellar 
populations. The dynamical properties of these galaxies, including the rotation 
speeds and evidence for high gas fractions and dark matter domination within 
the luminous extent of the galaxy, are consistent with those of previously studied 
faint galaxies in nearby groups. Further studies of the SDSS sample will allow 
characterization of low luminosity galaxies over the full range of environments in 
which they reside. 

Subject headings: galaxies: photometry, kinematics and dynamics 
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1. Introduction 

The Sloan Digital Sky Survey (York et al. 2000) galaxy redshift survey has an unprece- 
dented combination of large area, depth, and photometric quality, thanks to the combination 
of a large format camera (Gunn et al. 1998), high throughput multi-object spectrographs 
(A. Uomoto et al., in preparation), careful calibration procedures (Fukugita et al. 1996; 
Hogg et al. 2001; Smith et al. 2002), and an efficient series of data reduction and targeting 
pipelines (Lupton et al. 2001; Stoughton et al. 2002; Strauss et al. 2002; Blanton et al. 2003a; 
Pier et al. 2003). Recently Blanton et al. (2004a, hereafter B04) have searched the SDSS 
Second Data Release (DR2; Abazajian et al. 2004a) to identify a population of extremely 
low luminosity field galaxies. With absolute magnitudes M r ~ —13 to —16 (for the Hubble 
parameter h = ifo/lOOkms -1 Mpc -1 = 0.7 adopted throughout this paper), this population 
represents a range of luminosities that has previously been accessible to systematic study 
only in the Local Group and in nearby groups and clusters. 1 We have obtained Ha rotation 
curves of two of these objects, using the same methods that we are using for a study of the 
Tully-Fisher (1977, hereafter TF) relation of more luminous SDSS galaxies (J. Pizagno et 
al., in preparation, hereafter P04). Both objects appear to be "field" galaxies rather than 
satellites of brighter systems. Although the sample is small and the data quality limited, 
these measurements provide one of the first insights into the dynamical properties of this 
new population. 

2. Observations 

The galaxies in this paper were selected from the B04 sample according to their avail- 
ability during an observing run for the P04 project. As detailed in B04 and Blanton et al. 
(2004b), this low luminosity galaxy sample has been checked for contamination by double 
stars, errors in the photometric pipeline's automated deblending, and other complications 
that are unimportant for the great majority of SDSS galaxies but have a significant impact 
for low luminosity and low redshift systems. The two galaxies selected for follow-up obser- 
vations were J123654.9+013654.2 and J091858.60+581407.7. For brevity, we will hereafter 
refer to these as Galaxy 1 and Galaxy 2, respectively. 

Distances for galaxies in the B04 sample are estimated from their redshifts assuming 
Hubble flow and a model for the local peculiar field based on the IRAS 1.2-Jy redshift survey 



1 The field galaxy survey of Schombcrt, Pildis, & Edcr (1997), based on visual inspection of POSS-II 
plates and HI confirmation from Arecibo, includes some systems in this luminosity range. 
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(Willick et al. 1997). Details of this procedure are described by Blanton et al. (2004b). 
Because these galaxies have low redshifts (heliocentric values of 594 kms" 1 and 1157kms~ 1 , 
respectively), peculiar velocity uncertainties dominate the distance uncertainties, which in 
turn dominate the uncertainties in luminosity and dynamically inferred mass-to-light ratios in 
our discussions below. The Blanton et al. (2004b) procedure yields distances of 9.9 ±2.9 Mpc 
and 20.0 ± 2.2 Mpc to Galaxies 1 and 2, for h = 0.7. The la error bars are derived from the 
peculiar velocity probability function, which includes a 150 kms -1 local dispersion around the 
smooth velocity field predicted from the IRAS 1.2-Jy galaxies. The corresponding Petrosian 
r-band absolute magnitudes are M r = — 13.9+oJ and M r = — 14. 7^2 (again for h = 0.7). 
Neither galaxy has a bright neighbor (M r < —19.25) within 0.7 Mpc projected separation 
and 1000 km s -1 redshift separation. Galaxy 1 is close to NGC 4536 in projection, but the 
redshift difference, 594 km s -1 vs. 1808 km s -1 (Grogin, Geller, & Huchra 1998), makes a 
physical association unlikely. We conclude that both of these objects are "field" galaxies 
rather than satellites of brighter parents. 

Table 1 lists the magnitudes, redshifts, and distances for Galaxies 1 and 2 along with 
the g — r colors and axis ratios, exponential disk scale lengths, and rotation speed parameters 
described below. The top panels of Figure 1 show i-band images from the SDSS. A search 
through the NASA Extragalactic Database (NED) shows an HI observation (from Matthews 
& van Driel 2000) for the galaxy FGC 1475, which can be identified as our Galaxy 1 on the 
basis of angular proximity (12" difference between the NED and SDSS positions), similar 
axis ratios, and similarity between the Matthews & van Driel (2000) and SDSS redshifts. 

Spectroscopic observations were made with the CCDS long-slit spectrograph, set up to 
observe redshifted Ha, at the 2.4 meter MDM telescope during the night of 13 April 2004. 
We used a 2" slit width with a 600 lines/mm grating in the second order producing a 0.41 
A/pixel dispersion, and a 0.41 "/pixel spatial scale. The seeing during the night was 1.4 — 2". 

The Ha emission lines typically had a total signal-to-noise ratio of 6 — 15, with an 
intensity weighted velocity centroid uncertainty of ~ 8 km/sec. The dispersion axis was 
aligned to be perpendicular to the columns of the CCD to an accuracy of 0.1 A, as judged 
from the telluric lines. This alignment ensures that the Ha recession velocity is measured to 
the accuracy of the intensity weighted velocity centroids. The details of the data reduction 
and rotation curve extraction can be found in P04. 

The bottom panels of Figure 1 show the linearized and flat-fielded spectra of both 
galaxies. Figure 2 presents the extracted rotation curve measurements as observed veloc- 
ity relative to the continuum center vs. position along the slit in arc-seconds. Following 
Courteau (1997) and P04, we fit the rotation curve measurements with an arc-tangent func- 
tion, which has a minimal number of free parameters while still describing the global shape 
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of typical galaxy rotation curves quite well. Specifically, we use a Levenburg-Marquardt x 2 
minimization routine (Press et. al. 1992) to fit the data with the functional form 



where Vq is the central velocity, V c i IC is the asymptotic circular velocity, r is the position 
along the slit, r is the center of the rotation curve (where V = Vb), and r t is the turnover 
radius at which the rotation curve begins to flatten. The parameter uncertainties are derived 
from the covariance matrix returned by the Levenburg-Marquardt routine. For Galaxy 1, 
we allow all four parameters (Vb, V C i rc , r , r t ) to vary. For Galaxy 2, the best-fit arc-tangent 
parameters have very large V C \ TC and r t , so that they effectively describe a straight line over 
the region covered by the data points. We therefore fit a straight line to the rotation curve 
instead of an arc-tangent function, simplifying the determination of best-fit parameters and 
uncertainties. Since Galaxy 2's rotation curve is more extended on one side, we fix the line's 
intercept such that the line goes through the location of the conintuum center at the center 
of mass velocity. 

Figure 2 shows the arc-tangent and linear fits to the data points as smooth solid curves. 
For Galaxy 1, the outermost data points probe the turnover region of the rotation curve, while 
for Galaxy 2 they are still linearly rising. Clearly these data do not yield good constraints 
on the asymptotic circular velocity of a flat rotation curve. However, we find in P04 that 
the velocity V^d defined by the value of the arc-tangent fit at the radius of the outermost 
data point often provides a useful measurement of rotation speed even when the observed 
rotation curve is still rising at this radius. In particular, we find that if we use V cn d as the 
measure of rotation speed in the TF relation, then the observed TF scatter does not increase 
significantly when we include galaxies that have rising rotation curves in the sample. While 
the value of V^d is a lower limit to the asymptotic circular velocity, it can be much more 
robustly measured than V C i rc from truncated rotation curves like those in Figure 2. 

We correct V en d for inclination by using the GALFIT program (Peng et. al. 2002) to fit 
inclined exponential disks to the i-band galaxy images shown in Figure 1. The disk axis ratios 
and exponential scale lengths are listed in Table 1; formal uncertainties are ~ 0.01 in b/a 
and ~ 0.3" in -R e x P , though the systematic errors associated with assuming an exponential 
disk model are probably larger. The inclination corrected velocity is 



where 0.19 is the assumed intrinsic axis ratio for an edge-on disk and b/a is the measured 
i-band axis ratio. Different ranges for the intrinsic axis ratio vary by 0.10 to 0.20 depending 




(1) 




(2) 
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on galaxy type (Haynes & Giovanelli 1984). We chose 0.19, typical for spiral galaxies, and 
note that the range in intrinsic axis ratios causes a small variation (typically ~lkms _1 ) in 
Knd- GALFIT yields a position angle for Galaxy 2 that differs from the SDSS DR2 value 
by ten degrees, consistent with the visual impression of slight misalignment in Figure 1. We 
therefore apply a small (~ 1 kms" 1 ) correction to the measured V en d assuming a tilted-ring 
model as described by Beauvais & Bothun (2001; see their eq. 3). Table 1 lists inclination- 
corrected values V* nd including this slit misalignment correction for Galaxy 2. 

Matthews & van Driel (2000) report an HI velocity width for Galaxy 1 of = 
94 kms -1 , where the subscript c indicates a correction for instrumental resolution. They do 
not give an error bar on though the total signal-to-noise ratio of the HI line is 4.8, 

so the uncertainty in the width is probably not negligible. Kannappan, Fabricant, & Franx 
(2002) compare HI linewidths to optically measured, maximum rotation speeds and report 
a correlation 

W 50 = 19(±6) + 0.90(±0.03)(2\/ max ) , (3) 

in kms -1 units. Equation (3) predicts V max = 42 ± 5 kms -1 , compared to our measured 
value of V cnd = 34.8 ± 3.8 kms" 1 . (Note that the inclination correction is negligible for this 
galaxy.) The two measurements thus appear to be consistent within the observational errors 
(allowing for a few percent uncertainty in the measured W^c), though the somewhat higher 
value inferred from the HI data could reflect a continuing rise of the rotation curve beyond 
the radius probed by our Ha measurements. 



3. Discussion 

How do our measured rotation speeds compare to expectations based on the TF relation 
defined by more luminous galaxies? In P04, we use a sample of 170 galaxies with —22 < 
M r < —18 to measure the forward (M vs. log V) and inverse (log V vs. M) TF relations in 
the SDSS bands. For the r-band inverse relation, we fit 

rj = a(M r - M rfi ) + b (4) 

with i] = log 10 V^ nd and the constant M r $ = —20.873 chosen to yield minimal correlation 
between the errors in a and b. We find a = -0.143 ± 0.018 and b = 2.202 ± 0.004, yield- 
ing predicted values of V^ nd for Galaxies 1 and 2 of V pre & = lG.Olj^kms -1 and Vp re d = 
20-9^12 kms -1 , respectively. The uncertainties in M r and in the TF slope both contribute 
significantly to the errors on Vp rc d, and both contributions are non-linear. We have com- 
puted upper la - error bars by separately varying M r by —la and a by +la and adding 
the corresponding values of SV prc d in quadrature; we follow an analogous procedure for the 
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lower la error bar. The predicted rotation speeds are well below the measured values of 
34.8 ± 3.8 kms^ 1 and 30.9 ± 7.2 kms -1 , though marginally compatible in the latter case be- 
cause of the large observational errors. Dotted curves in Figure 2 show rotation curves with 
the same parameters as the solid curves but scaled by {V pre ,i/V^ nd ). Despite the uncertainties 
in the measurements, it is clear that the galaxies are rotating substantially faster than pre- 
dicted by an extrapolation of the inverse TF relation to these faint magnitudes. To match 
the observed rotation speeds while maintaining the P04 normalization at M r ss — 21 would 
require an inverse TF slope of a ks —0.1, compared to the slope of —0.143 measured in the 
bright galaxy regime. 

We can use our measurements to estimate dynamical mass-to-light ratios within the 
radius of the outermost Ha data point in each galaxy. We use the simple mass estimator 

M = ( ^ nd)2fiend , (5) 
G 

ignoring possible corrections for asymmetric drift or non-circular motions, which would be 
small compared to our observational error bars. The values of i? en d are 19 and 11 arc-seconds 
for Galaxies 1 and 2, respectively, corresponding to 1.7 and 1.8 times the exponential scale 
lengths measured by GALFIT. At our estimated distances of 9.9±2.9 Mpc and 20±2.2 Mpc, 
the physical values of i? e nd are 0.91 ± 0.27 kpc (Galaxy 1) and 1.07 ± 0.12 kpc (Galaxy 2). 
We measure the fraction of r-band light within R cn d for each galaxy and correct the values 
of M r in Table 1 accordingly, by 0.24 and 0.14 magnitudes. We do not apply any internal 
extinction corrections. Adopting M r:Q = 4.67, we find M/L = 12. M Q /L Q for Galaxy 1 
and 4.8l|i M e /L Q for Galaxy 2. Because M/L oc V^ nA /d and the fractional uncertainties in 
these quantities are substantial, we compute the upper la error bar by separately varying 
d by — la and V cn d by +la and adding the two changes 5M in quadrature, and we follow 
an analogous procedure for the lower la error bar. The M/L value for Galaxy 1 is much 
higher than the value ~ 1.3M /L expected for a stellar population with the observed g — r 
color of the galaxy (Bell et al. 2003), so this system is strongly dominated by dark matter 
within -Rend- The case for dark matter domination in Galaxy 2 is less clear because of the 
large uncertainty in M/L, though the central value again corresponds to a large ratio of dark 
matter to stellar mass within -R en d- (The Bell et al. [2003] models imply a stellar mass-to- 
light ratio of ~ O.9M /L for Galaxy 2's g — r color.) While all disk galaxies become dark 
matter dominated at sufficiently large radii, these low luminosity systems appear to have 
high dark matter fractions even within two disk scale lengths. 

To put our measurements in the context of previous results, Figure 3 plots inclination 
corrected circular velocities against -B-band absolute magnitude for our galaxies (filled circles 
with error bars) and for galaxies with HI line widths and Mb > — 15 from Cote, Carignan, 
& Freeman (2000; open squares), Begum & Chengalur (2004; open circles), and HI line 



-7- 



widths from Carignan & Freeman (1988; open triangle). For our galaxies, we use V* nd as 
the indicator of circular velocity, and we convert the (AB) g-band luminosity to the (Vega) 
.B-band luminosity using the galaxy g — r color and M B = M g + 0.365 + 0.46 [(p — r) — 0.78], 
obtained using the i^-correction code of Blanton et al. (2003b). For Cote, Carignan, & Free- 
man (2000) we use the value of V c listed in their tables, and for Begum & Chengalur (2004) 
and Carignan & Freeman (1988) we estimate V c visually from the plotted rotation curves. 
We have not attempted to put observational error bars on the literature data points because 
this would require a detailed assessment of the uncertainties in the distance measurements, 
which use a somewhat different method in each case. Cote, Carignan, & Freeman (2000) 
observe galaxies in the Sculptor and Centaurus A groups and assign group distances based 
on a variety of indicators (see Cote et al. 1997 and references therein). Begum & Chengalur 
(2004) assign both of their galaxies to the NGC 4696 group on the basis of proximity and 
adopt a distance based on the brightest stars (Huchtmeier et al. 2000). Carignan & Freeman 
(1988) use a combination of group assignment and brightest stars to infer the distance of 
DDO 154. 

The solid line in Figure 3 shows the 5-band inverse TF relation extrapolated to the low 
luminosity regime. We convert the g-band relation of P04 to .B-band using Mb = M s + 0.19, 
appropriate for a galaxy with g—r = 0.4. Dotted lines show relations with the slope varied by 
±lo". For bright galaxies, the intrinsic scatter about the mean relation is 0.07 dex in log 10 V^ nd 
(P04), smaller than the 0.1 dex observational error bar on Galaxy I. As noted earlier (on 
the basis of r-band data), Galaxy 1 is rotating substantially faster than the extrapolated TF 
relation predicts, while the large error bar on V cn d leaves Galaxy 2 marginally consistent with 
the TF extrapolation. The rotation speed of Galaxy 1 appears reasonably in line with that 
of the other dwarf galaxies, while the rotation speed of Galaxy 2 is noticeably low, perhaps 
because the optical rotation curve is still rising steadily at our outermost data point. All of 
these systems are rotating faster than the TF extrapolation predicts; equivalently, they are 
fainter than predicted given their rotation speeds. 

McGaugh et al. (2000) show that low luminosity galaxies are fainter than predicted 
by the (forward) TF relation defined by bright galaxies, consistent with the result shown in 
Figure 3. However, faint galaxies are often gas rich, and they show that adding the gas masses 
inferred from HI observations to the stellar masses leads to a linear relation between log M bar 
and log 14 over four decades in baryonic mass M bar , extending down to a few xlO 7 M . For 
our Galaxy 1, Matthews & van Driel (2000) report an HI flux integral of 1.24 Jykms -1 , 
which implies an HI mass to optical luminosity ratio of ~ 0.95M Q /L Q in r-band, or roughly 
equal mass in neutral gas and stars. This ratio is below the median of ~ 2M /L found (in 
V-band) by Pildis, Schombert, & Eder (1997), but within the range spanned by their data. 
Simply doubling the luminosity of Galaxy 1 and thus shifting it rightwards by 0.75 mag in 
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Figure 3 would not move it onto the P04 inverse TF relation. However, McGaugh et al. 
(2000) find a steep slope for the (forward) baryonic mass TF relation, and the baryonic mass 
of 7.5 x 1O 7 M inferred for Galaxy 1 at a distance of 9.9 Mpc agrees respectably with the 
value M bar = 3O.5(/i/O.7)- 2 (\/ c /kms- 1 ) 4 M = 5.3 x 1O 7 M predicted by their fitted relation 
for V c = 35kms _1 , given the substantial error bars on both the predicted and measured 
values. 

The luminosity-velocity relations and mass-to-light ratios of low luminosity galaxies 
offer clues to the role of supernovae and photoionization in regulating star formation, since 
these feedback effects are generally expected to be stronger in lower mass halos (see, e.g., 
Dekel & Silk 1986; Quinn, Katz, k Efstathiou 1996; Thoul & Weinberg 1996; Benson et 
al. 2002; Somerville 2002). The influence of feedback may be different in galaxies that are 
central objects of their parent dark matter halos and in galaxies that are satellites in more 
massive halos — semi-analytic models and hydrodynamic simulations predict systematically 
different properties for these two populations (Kauffmann, White, & Guiderdoni 1993; Cole 
et al. 1994; Somerville & Primack 1999; Berlind et al. 2003; Zheng et al. 2004). Most of the 
very low luminosity galaxies studied to date have been discovered or investigated because 
they are members of groups or clusters defined by brighter galaxies. The two galaxies studied 
in this paper are not satellites of bright neighbors, but their properties are roughly in line with 
those of satellite dwarfs. The B04 sample drawn from the SDSS provides the opportunity to 
study such faint galaxies in detail over the full range of environments in which they appear. 
Our present investigation represents a start on this broader program. 
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Table 1. Galaxy Properties 



Galaxy 


r 


g-r 


cz 


d ( Mpc) 


M r 


-Rexp (") 


b/a 


Vend 


^prcd 


J123654.9+013654.2 
J091858.60+581407.7 


16.11 
16.84 


0.41 
0.31 


594 
1157 


9.9 ±2.9 
20.0 ±2.2 


-13.9121 
-14.7t°;2 


11.2 

6.2 


0.20 
0.40 


34.8 ±3.8 

30.9 ± 7.2 


i6.ot 6 5 i 

20.9tti 



Note. — Velocity units (for cz, V en( j, and V pre d) are km s" 1 . Distances and absolute magnitudes are computed assuming 
Hq = 70kms~ 1 Mpc -1 , using a model of the local peculiar velocity field as described in the text. All the magnitudes are 
corrected for Galactic foreground extinction. 



Fig. I. — Images and spectra of Galaxy 1 (left) and Galaxy 2 (right). Upper panels show 
the SDSS i-band images with North down and East to the left. The small vertical line has 
a length of 20 arc-seconds, and the long line through each galaxy shows the slit position 
angle. Lower panels show the flat-fielded, linearized spectra from MDM. These panels are 
40 arc-seconds high and 61 A wide, centered at 6582 A. 
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Fig. 2. — Rotation curves of Galaxy 1 (left) and Galaxy 2 (right). Points with error bars 
show the flux-weighted velocity centroid of the Ha line as a function of position along the 
slit. For Galaxy 1, the solid curve shows the arc-tangent fit to the data (eq. 1). For Galaxy 
2, we use a straight fit. Our measure of circular velocity is the value V en d of the smooth 
fit at the radius of the outermost data point. Dotted curves show the arc-tangent or linear 
fits scaled by Vp rc d/Knd) where V^ re( j is the value predicted by extrapolating the inverse TF 
relation of P04 to the absolute magnitudes of these galaxies. Note that these plots are not 
corrected for inclination. 
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Fig. 3. — Inclination corrected circular velocity vs. 5-band absolute magnitude for our 
galaxies (filled circles with error bars) and for low luminosity galaxies with HI data taken from 
Cote, Carignan, & Freeman (2000, open squares), Begum & Chengalur (2004, open circles), 
and Carignan & Freeman (1988, open triangle). The solid line shows the extrapolation of 
the 5-band inverse TF relation to the faint galaxy regime. The 5-band inverse TF relation 
is found by converting the g-band inverse TF relation of P04. Dotted lines show this relation 
with the slope changed by ±la. 



